Statistical and dynamical fluctuations of Binder ratios in heavy ion collisions 
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Higher moments of net-proton Binder ratio, which is suggested to be a good experimental mea- 
surement to locate the QCD critical point, is measured in relativistic heavy ion collisions. We 
firstly estimate the effect of statistical fluctuations of the third and forth order Binder ratios. Then 
the dynamical Binder ratio is proposed and investigated in both transport and statistical models. 
The energy dependence of dynamical Binder ratios with different system sizes at RHIC beam scan 
energies are presented and discussed. 

PACS numbers: 25.75.Gz,25.75.Nq 



I. INTRODUCTION 



One of the main goals of current relativistic heavy ion 
collisions is to map the QCD phase diagram At 
vanishing baryon chemical potential /is = 0, finite tem- 
perature Lattice QCD calculations predict that a cross- 
over transition from hadronic phase to the Quark Gluon 
Plasma (QGP) phase will occur around a temperature 
of 170 - 190 MeV 0,11. QCD based model calculations 
indicate that the transition could be a first order at large 
fJ-B [J- The point where the first order phase transition 
ends is the so-called QCD Critical Point (QCP) [1, [Ij. 
Attempts are being made to locate the QCP both ex- 
perimentally and theoretically Q. Lattice QCD calcula- 
tions at finite hb face numerical challenges in comput- 
ing. Thus the location of the QCP are highly uncertain 
in theoretically. In experimental aspect, the RHIC beam 
energy scan program Q has been motivated to search 
for the QCP in experiment. By decreasing the collision 
energy down to a center of mass of 5 GeV, RHIC will be 
able to vary the baryon potential from hb ~ to 500 
MeV. 

Fluctuations of conserved charges, which behave dif- 
ferently between the hadronic and QGP phase, are gen- 
erally considered to be sensitive indicators for the tran- 
sition d, [l^l- The singularity at the QCP, at which the 
transition is believed to be second order, may cause en- 
hancement of fluctuations if fireballs created by heavy 
ion collisions pass near the critical point during the time 
evolution [ll|. It has been shown that near the critical 
point, the density-density correlator of baryon number 
follow the same power law behavior as the correlator of 
the sigma field which is associated with the chiral or- 
der parameter (ill . [T2j . Therefore, the baryon number is 
considered as an equivalent order parameter of formed 
system in nuclear collisions. In experiment, net-proton 
multiplicity distribution is much easier to measure than 
the net-baryon numbers. Theoretical calculations have 
shown that in QCD with exact isospin invariance, the 
relevant corrections due to isospin breaking are small 
and the net-proton fiuctuations can reflect the singular- 
ity of the baryon number susceptibility as expected at the 
QCP [H]. Hence, the net-proton number is used as the 



order parameter in current heavy ion experiment |i4| . 

It is suggested recently that the Binder-like ratios 
are good identification of critical behavior in relativis- 
tic heavy ion collisions [T5l - [T7| . In statistical physics, the 
Binder ratios arc defined as normalized raw moments of 
order parameter. The third and forth order Binder ratios 
can be written as 
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where M is the order parameter. The Binder ratio is 
a new observable in heavy ion collisions. The differ- 
ence between Binder ratios and the well-known higher 
moments ITsIl is that Binder ratios are the normalized 
higher raw moments of the order parameter while the 
higher moments use central moments. 

The universality argument indicates that the static 
critical exponents of the second order phase transition are 
determined by the dimensionality and symmetry of the 
system. The QCD critical point of deconfinement phase 
transition belongs to the same universality class as liquid- 
gas phase transition and the 3D-Ising model [ll|, [l^, [l^ . 
Its universal critical properties are discussed to be valid 
in various of models and relevant to heavy ion colli- 
sions [ill, 113, im , in particular the event-by-event fluc- 



tuations of baryon numbers [11 



In the calculations of the 3D-Ising model[l5|. Binder 
ratios of B3 and B4 as a function of temperature show 
a step jump from a lower platform to a higher one near 
the vicinity of critical point by choosing the magnet of 
spin as the order parameter. Therefore, if the formed 
system in heavy ion collisions reaches the critical point 
and the freeze-out curve is close to the transition line, the 
step function liked behavior of Binder ratios in the Ising 
model may serve as a good probe of QCP in current heavy 
ion collisions, where critical incident energy is difficult to 
assign precisely in priori. 

In the mean time, the effect of trivial statistical fluctu- 
ations (22I I due to insufficient number of particles should 
be studied and properly eliminated in higher moment 
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calculations. Therefore, we should discuss the statisti- 
cal contributions from the measured fluctuations firstly, 
then we could identify the dynamical part which is more 
relevant to the critical point of the QCD phase transition. 

In this paper, we firstly investigate the statistical and 
dynamical fluctuations of net-proton Binder ratios by us- 
ing the AMPT and THERMINATOR models. Then, the 
energy dependence of dynamical Binder ratios in Au -I- 
Au collisions at various RHIC beam scan energies are 
studied. 

In our analysis, two versions of a multi-phase trans- 
port (AMPT) model [H are used. One is the AMPT 
default and the other one is the AMPT with string melt- 
ing. In both versions, the initial conditions arc obtained 
from the heavy ion jet interaction generator (HIJING) 
model, and then the scattering among partons is given 
by the Zhangs parton cascade (ZPC) model. In the 
AMPT default model, the partons recombine with their 
parent strings when they stop interacting, and the re- 
sulting strings are converted to hadrons using the Lund 
string fragmentation model, whereas in the AMPT model 
with string melting, quark coalescence is used in combin- 
ing partons into hadrons. The dynamics of the hadronic 
matter is described by the ART model. The THERMI- 
NATOR statistical model [11] is a Monte Carlo event 
generator designed for studying of particle production in 
relativistic heavy ion collisions from SPS to LHC ener- 
gies. It implements thermal models of particle produc- 
tion with single freeze out. 

In order to make our calculations convenient for com- 
parison with the RHIC beam energy scan data, we choose 
the midrapidity {\y\ < 0.5) region with transverse mo- 
mentum 0.4 < pT < 0.8 GeV/c. This phase space is 
where the STAR experiment can do the particle identi- 
fication for proton numbers with its main tracking de- 
tector - the Time Projection Chamber [T!] . The number 
of events used in this analysis is around 6 million. This 
statistics is needed for the calculation of the dynamical 
Binder ratios of net-proton to ensure the statistical errors 
under control. 
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FIG. 1: Binder ratios of (a) and B4 (b) as a function 
of number of participants from AMPT default (solid cir- 
cle), AMPT with string melting (open circle), and THER- 
MINATOR (open square) models in Au + Au collisions at 
y/SNN = 200 GeV. The dashed lines represent the correspond- 
ing statistical fluctuations. 
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where A = (iVp) — (Np) is the average number of net- 
protons, and /i = {{Np) + (Np)) /2 is the mean value 
of protons and antiprotons in the event sample. More 
details of the calculation of this formula could be found 
in the appendix. 

In Fig. 1 (a) and (b), we show the results of and 
Bi as a function of number of participants (iVpart) from 
AMPT default (solid circle), AMPT with string melt- 
ing (open circle), and THERMINATOR (open square) 
models in Au -I- Au coUisions at ^/sNN = 200 GeV, re- 
spectively. For comparison, the statistical fiuctuations of 
the Binder ratios, calculated from Eq. (2), are presented 
as dashed lines. We can see that in both transport and 
statistical models, the statistical fiuctuations give main 
contributions to the Binder ratios. It shows that the 
influence of statistical fluctuations are not negligible in 
the measurement of net-proton Binder ratios at RHIC 
energy. 



II. STATISTICAL FLUCTUATIONS OF BINDER 
RATIOS 

In the measurement of the net-proton Binder ratios, 
finite number of protons and antiprotons will cause non- 
negligible statistical fiuctuations. If the produced pro- 
tons and antiprotons are two independent Poisson-like 
distributions (2^ . [25| , the net-protons then obey a Skel- 
1am (SK) distribution 

According to the definition of Eq. (1), the statistical 
fiuctuations of the net-proton Binder ratios can be di- 
rectly deduced from the Skellam distribution 



III. DYNAMICAL NET-PROTON BINDER 
RATIOS 

In the THERMINATOR model, it is well-known that 
the fiuctuations are thermal. From Fig. 1, we observe 
it gives a good agreement with the Skellam statistical 
fiuctuations. It is difficult to disentangle purely statis- 
tical effects from thermal fluctuations which follow the 
physics of a hadron resonance gas. Since neither of them 
is associate with the QCP behavior, we suggest to elimi- 
nate these statistical or thermal fiuctuations in order to 
get the dynamical part. 

As shown in section II, the statistical fiuctuations of 
Binder ratio can be expressed by Eq. (2) given proton 
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FIG. 2: The third (a) and forth (b) order dynamical Binder 
ratios as a function of A'^part from AMPT default (solid cir- 
cle), AMPT with string melting (open circle), and THER- 
MINATOR (open square) models in Au + Au collisions at 
y^sWN = 200 GeV. 



and antiproton obey independent Poisson distributions. 
We define the so-called dynamical Binder ratios as, 



m 0.1 




Bs^dyn — — B^^stat 

Bi^dyn = B4 — Bi^stat- (3) 

We suggest to measure these dynamical Binder ratios 
instead of the original definition given by Eq. (1) in rcl- 
ativistic heavy ion experiment. 

The dynamical net-proton and as a function 
of iVpart from AMPT default (solid circle), AMPT with 
string melting (open circle), and THERMINATOR (open 
square) in Au -I- Au collisions at ^Jsnn — 200 GcV arc 
shown in Fig .2 (a) and (b), respectively. We find that 
both the third and forth order Binder ratios from THER- 
MINATOR are zero at all ccntralities. This is because 
that THERMINATOR is based on the hadron resonance 
gas model and the produced net-protons in the final state 
obey the Skellam distribution [231 ■ While, in transport 
models, both dynamical B^ and B^ are larger than zero in 
peripheral collisions, then tend to be zero in central col- 
lisions. The results from AMPT string melting are larger 
than that from the default model. This is due to different 
mechanisms of hadronization scheme used for finite state 
particles in different versions of AMPT models. 



IV. ENERGY DEPENDENCE OF DYNAMICAL 
BINDER RATIOS IN TRANSPORT MODELS 

Based on the Wolf algorithm the 3D-Ising model 
calculations show that the third and forth order Binder 
ratios exhibit a step function behavior near the critical 
temperature [l^. In experiment, if the freeze-out curve 
is close to the QCD critical point along the transition 
line and the critical fluctuations are survived in the final 
state, one may expect to observe this behavior exper- 
imentally. For RHIC beam energy scan program, if we 
scan all the existing energies and observe two platforms of 
Binder ratios at low and high energy regions, respectively, 



FIG. 3: Energy dependence of the dynamical Bi (left panel) 
and Bi (right panel) at RHIC energies from AMPT default 
model (upper panel) and string melting (lower panel), respec- 
tively. The nine different symbols represent nine collision 
sizes, which goes from most peripheral (70-80% central) to 
most central (0-5% central) collisions. The lines are used only 
to guide eyes. 

the critical point is most probably in the intersection of 
different size systems between the two platforms. 

In the upper panels of Fig. 3 (a) and (b), we show 
the energy dependence of the dynamical i?3 and -B4 at 
six RHIC energies, 7.7, 9.2, 11.5, 39, 62.4, and 200 GeV 
from AMPT default model. The nine different symbols 
represent nine collision sizes (denoted by centralities in 
experiments). From the figure, we can sec that, in all 
beam energies, when centrality goes from most peripheral 
(70-80% central) to most central (0-5% central) collisions, 
the dynamical Binder ratios decrease and are close to 
zero in the most central collisions. It means that both 
dynamical iJs and B/^ are system size dependent. In the 
calculations of 3D-Ising model [l^ , the Binder ratios of 
different sizes should intersect to each other, i.e. there is 
a fix point for the system near the critical point. From 
7.7 GcV to 200 GeV, wc observe no fixed point from the 
AMPT default model. The lower panels of (c) and (d) 
for the string melting version give similar results. 

Therefore, there is no step function behavior observed 
in both two versions of the AMPT models. This is un- 
derstandable since there is no QCD critical mechanism 
implemented in these transport models. 

V. SUMMARY AND OUTLOOK 

In this paper, the statistical and dynamical Binder ra- 
tios of net-proton are studies in Au -|- Au collisions at 
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RHIC energies. Using transport and statistical models, 
it is shown that statistical fluctuations are not negligible 
in the measurement of higher Binder ratios in relativistic 
heavy ion collisions. 

In order to obtain a clean signature which may be re- 
lated to the QCD critical point, we suggest to use the dy- 
namical Binder ratio in experimental measurement. The 
dynamical net-proton Binder ratio is found to be zero in 
the THERMINATOR model but larger than zero in pe- 
ripheral collisions in AMPT model. The energy depen- 
dence of dynamical Binder ratios with different system 
sizes shows no step function behavior either in AMPT 
default or string melting models, which is consistent with 
expectation that no QCD critical point physics is imple- 
mented in these models. 

It is interesting to sec whether the step function behav- 
ior of the dynamical Binder ratios could show up in the 
coming high energy collisions at RHIC, SPS, and FAIR 
experiments, where the critical incident energy of QCD 
phase transition may be covered. Our model study can 
serve as a background study of the behavior expected 
from known physics effects for the experimental search 
for the QCD critical point. 
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VII. APPENDIX: BINDER RATIOS FROM 
SKELLAM DISTRIBUTION 

Given the distributions of proton and antiproton are 
independent Poisson distributions with mean values are 



(Np) and (Np), the net-proton will follow a Skellam dis- 
tribution. If we define the net-proton as M = Np — Np, 
then the probability distribution function of M is 

f{M; {Np),{Np}) 

where /|m| ^2-^ (Np) (Np)^ is the modified Bcssel func- 
tion of the first kind. 

The 71*'' moment of M, which is defined as (M") = 
IZo M'^fiM; {Np), {Np))dM, can be calculated from the 
above distribution function. We obtain 

(Af) = A 

(M^) = ^? + 6/iA + A 

(M'') = A" -f 12/iA2 4A2 12Ai2 2/^, 
where A = {Np) — {Np) is the average number of net- 
protons, and = {{Np) + {Np)) /2 is the mean value of 
proton and antiproton. 

By the definitions of the third and forth Binder ra- 
tios of Eq. (1), wc get the Binder ratios of the Skellam 
statistical distribution as 



_ < AP > _ + 6^iA + A 

- < M2 >3/2 " (A2+2//)3/2 ' 
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